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Abstract 

The estimate of the energy deposition rate (EDR) for neutrino pair annihilation has been carried 
out. The EDR for the neutrinos coming from the equatorial plane of a rotating neutron star is 
calculated along the rotation axis using the Cook-Shapiro- Teukolsky (CST) metric. The neutrino 
trajectories and hence the neutrino emitted from the disk is affected by the redshift due to disk 
rotation and gravitation. The EDR is very sensitive to the value of the temperature and its 
variation along the disk. The rotation of the star has a negative effect on the EDR, it decreases 
with increase in rotational velocity. 
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I. INTRODUCTION 



Gamma Ray Bursts (GRB), first discovered in the late 1960s by US military satellites 
and its possible connection with neutrino production in compact stars is a field of high 
current interest. The central engine of GRBs is still shrouded in a mist, although several 
models has been proposed. GRBs are separated in two classes; long duration bursts (long 
GRB) which lasts from 2 second to several minutes, with average duration of 30 seconds 
and short duration bursts (short GRB), with burst duration from few milliseconds to 2 
second, average being 0.3 seconds. Of late, the detection of afterglows in short GRBs and 
also precise localization of different short GRBs have provided new inputs for the study of 
these phenomena. The light curves ^ and polarization j^] of gamma ray burst afterglows 
indicate anisotropic central engines. These observations have been considered as evidence 
that GRBs are highly beamed. 

Neutron stars are objects formed in the aftermath of supernova. The central density of 
these stars can be as high as 10 times that of normal nuclear matter. At such high density, 
any small perturbation, e.g. spin down of the star, may trigger a phase transition from 
nuclear to quark matter system. As a result, the neutron star may convert to a quark 
star or a hybrid star with a quark core j^, O, Q]. It has been shown [71] that such a phase 
transition 8| produces a large amount of high energy neutrinos. These neutrinos (and 
antineutrinos) could annihilate and give rise to electron-positron pairs through the reaction 
vv — > e + e~ . These e + e~ pairs may further give rise to gamma rays which may provide a 
possible explanation of the observed GRB. Furthermore, the rotating neutron stars have 
been shown |9[ to produce the observed beaming effect. At present, it is necessary to have a 
better understanding of the energy deposition in the neutrino annihilation to e + e~ process 
in the realistic neutron star environment. 



Motivated by the delayed explosion of Type II supernovae, the 



1DR due to the reaction 



vv — > e + e~ in the vicinity of a neutron star have been calculated [1PJ,1H| based on Newtonian 
gravity, i.e. (2GM/c 2 R) << 1. Goodman et al. Hj pointed out that the neutrino pair 
annihilation rate can be seriously altered by the gravitational effects. The effect of gravity 
was incorporated in refs. 0, Q], but only for a static star. Slow rotation on this calculation 



was introduced by Prasanna and Goswami 



14j. GR effect on energy deposition rate for 



neutrino pair annihilation near a black hole was studied in ref. 15|, lla] to explore the 
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possible engine for GRBs. Here the most probable candidate for the central engine was the 
accretion disk around a black hole. A detailed hydrodynamic simulation of the process was 
studied by Birkl et al. |17 ]. 

The relativistic effect consists of three factors: the gravitational redshift, the bending 
of the neutrino trajectories and redshift due to rotation. The EDR is enhanced by the ef- 
fect of neutrino bending, however the redshift due to disc rotation and gravitation reduces 
the EDR. Thus these effects further complicate the estimation of EDR. In order to discuss 
quantitatively the central engine of GRBs, we need a comprehensive study of the forma- 
tion, evolution, temperature dependence, geometry of the star, the mechanism of energy 
deposition etc., most of which are affected by the rotational and GR effect. 

In an earlier paper [l8| we have shown that the path of neutrino is important for the 
study of EDR near a massive object. We have done a complete GR calculation of the 
neutrino path for the most general metric describing a rotating star, and obtained its geodesic 
equation along the equatorial and polar plane. The minimum photosphere radius (MPR) 
has been calculated for various stars along these two planes. Our aim in this work is to 
study semianalytically the relativistic effect on the EDR above the equatorial plane for a 
star rotating along the polar axis. We calculate the EDR along the rotation axis for the 
neutrinos coming from the equatorial plane of the star. 

Our paper is arranged as follows. In the next section we fix our EOS and describe the 
geometry of the star. In section III we formulate the algorithm of the EDR calculation 
above the equatorial plane along the rotation axis. Next we present our results and finally 
we discuss the results and summarize our work. 

II. STAR STRUCTURE 

The structure of the star is described by the CST metric given by [3] 

ds 2 = -e (+p dt 2 + e 2a {dr 2 + r 2 d6 2 ) + e^- p r 2 sin 2 d(d<P - udt) 2 . (1) 

The four gravitational potentials, namely a,7,p and u are functions of 9 and r only. All 
the potentials have been solved for both static as well as rotating stars using the 'rns' code 



2d, 



21 



221 ] . For a rotating star, all the potentials become functions of both r and 9. 



Tabulated equations of state (EOS) are needed to compute the code numerically. In 
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this paper we have used hadronic EOS evaluated using the nonlinear Walecka model [23]. 
The Lagrangian in the model includes nucleons (neutrons and protons), electrons, isoscalar 
scalar, isoscalar vector and isovector vector mesons denoted by ipi, ip e , a, uo^ 1 and p a,fl , 
respectively. The Lagrangian also includes cubic and quartic self interaction terms of the a 
field. The parameters of the nonlinear Walecka model are meson-baryon coupling constants, 
meson masses and the coefficient of the cubic and quartic self interaction of the a mesons. 
The meson fields interact with the baryons through linear coupling. The uj and p meson 
masses have been chosen to be their physical masses. The rest of the parameters, namely, 
nucleon-meson coupling constants and the coefficients of cubic and quartic terms of the a 
meson self interaction are determined by fitting the nuclear matter saturation properties, 
namely, the binding energy/nucleon (-16 MeV), baryon density (po=0.17 fm~ 3 ), symmetry 
energy coefficient (32.5 MeV), Landau mass (0.83 m n ) and nuclear matter incompressibility 
(300 MeV). 

Using this tabulated EOS and a fixed central density, the rns code is solved to obtain all 
the gravitational potentials as a function of r and 9, which thereby defines the shape, mass, 
rotational velocity and other par ameters of the star. The shape of a fast rotating neutron 
star becomes oblate spheroid 19j. The star gets compressed along the z-axis, on the other 
hand along x and y-axes it bulges by equal amounts, the polar radius thus being smaller 
than equatorial radius. 

III. GR CALCULATION 

In this section we study the GR effect on the neutrino pair annihilation rate. The coor- 
dinate system is oriented in such a way that the equatorial plane lies along 9 = | and the 
polar plane along 9 = 0. We treat the equatorial plane as a disk from which the neutrinos 
are coming out and the EDR is calculated along the different radial points on the rotation 
axis. Here we are interested in the EDR via neutrino pair annihilation near the rotation 

It has been shown earlier [15J that in the absence of gravitation the 
9 dependence of the EDR is weak for small values of 9. Although the 9 dependence of the 
gravitational effect may not necessarily be small, here we calculate the EDR at 9 = and 
assume it to be approximately same for small values of 9. 



The energy deposited per unit volume per unit time is given as 10, 111] 
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where f u (ft?) is the number density of neutrino (antineutrino), v v is the neutrino velocity, 
and a is the rest frame cross section for the process vv —>■ e + e~ . The Lorentz invariant term 
is given by 

r F 



{a\V v — V v \e v e v ) 



DGi 2 



3tt 



with 



G F 2 = 5.29 x 10~ u cm 2 MeV- 2 ; 
D = l± Asin 2 9 w + 8sin 4 6 w . 



(3) 

(4) 
(5) 



sin6 w = 0.23. (+) sign is for electron neutrino-antineutrino pair and (— ) sign for other 
pairs. 

As the geometry of the equatorial plane is circular, we can decouple the energy and 
angular dependence. Thus the rate of energy deposition is 
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where, F(r), the angular integral, is given by 161 ] 
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We calculate the energy integral numerically, integrating it from to oo. Both the angular 
and energy integral has temporal component which are affected by the GR and rotational 
effect. T e ff is the effective temperature of the disc that is observed in the comoving frame. 
The value of T is the temperature of the disc that is observed at infinity. In the locally 
nonrotating frame a neutrino moves normally to the direction of the disk motion. In that 
case, the temperature suffers redshift both due to disk rotation and gravitation. Therefore 
the two temperature T and T e ff are related [la], as 

T eff 
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gtt, Qt<f> and are calculated from the metric. T is the Lorentz factor defined as T = ^_ v2 
with v being v — (A — u)rsinOe~ p , A the rotational velocity of the star. 



In an earlier paper [18j a detailed GR calculation of the neutrino path has been done. 
For the present coordinate system, the 4- momentum [16( is given by the equation, 

9^P U + Vm 2 = 0, (9) 

where /i m is the rest mass of the particle and p M = A being the afflne parameter. After 
a bit of algebra the final geodesic equation can be written as, 

2al A(r,6) l2 e^)r 2 sm 2 ^ 2 e 2 ^6 ]2 

L A(r,0) +cj e 2 ° 1 v 7 r 2 s%nW 

-e^ +p \l - ubf + 9 ^ 9/ , e^ +3p ) = 0. (10) 

This equation can be solved using the potentials obtained from rns code to obtain a mini- 
mum radius r = Rmpr, the minimum photosphere radius, below which a massless particle 
(neutrino) emitted tangentially to the stellar surface (0r = 0) would be gravitationally 
bound. This is also the minimum radius from which the neutrinos can come out from the 
disk. 

From an angle 6 V at a point (r, 0), we can trace back a trajectory of the neutrino to its 
emission point on the disc. Thus the temperature To depends on 9 V , and therefore T e ff 
appears in the integration of Q v in the angular integral. The neutrinos are coming out of 
the disk and depositing their energy along the rotation (symmetry) axis of the disk. The 
minimum radius from which the neutrinos can come out is given by MPR. The outermost 
point from which the neutrinos comes out is from the surface of the disc {%. e. the surface 
of the star along the equatorial plane). Figure 1. shows the layout of our problem, where 
neutrinos coming out from the disk and are depositing their energy along the rotation axis. 
The angle from which the neutrinos are coming can be calculated from the figure. The angle 
8 m is the angle subtended by the MPR on the rotation axis, which is given by 

e m = tan- 1 ^^. (11) 
r 

where Rmpr is the MPR on the disc. Similarly 9m is the angle subtended by the surface of 
the star on the rotation axis 

M = tan -i^L, (12) 



where Rs U r is the radius of the outer surface of the disc i. e. equatorial radius of the star. 

The EDR is now calculated by integrating q over the proper volume. We calculate the 
amount of energy deposited over r and within 9 in the close vicinity of the rotation axis 15]. 
The point on the rotation axis from where we start our calculation is r = Rmpr- in this 
problem we are considering only the thermal neutrinos and our lower limit of integration 
is the MPR. Beyond the MPR the density of the star is not as high as that of the central 
region. Beyond the MPR the thermal neutrinos are transparent to the star and therefore 
the rate of deposition energy reabsorbed is not so large and we ignore it. 



IV. RESULTS 

We choose the central density of the star to be 1 x 10 15 gm/cm 3 , for which the Keplerian 
velocity is 0.61 x 10 4 s -1 . For comparison we have also done the calculation for a slow star 
(0.4 x 10 4 s _1 ). For the keplerian star the equatorial radius is 16km, which serves as the 
outermost point of the disk. The innermost point from which the neutrinos can contribute 
in the EDR is the MPR, which is 4km from the centre. The same two points for the slowly 
rotating star are 12km and 4.7km respectively [Lsl ]. 

First we assume an isothermal disk, i.e. there is no temperature gradient (TG) along 
the disk. Figure 2 shows the variation of the EDR along the rotation axis of the star for 
three different temperature (lMeV, 5MeV and lOMeV). Initially the EDR increases with 
distance and it reaches a maximum value at a distance of 16km and after that it falls off 
gradually. The EDR is maximum for central temperature Tq = lOMeV and minimum 
for Tq = lMeV, which shows that the EDR increases with increase in temperature of the 
disk. For Tq = lOMeV the maximum EDR is at a point about 16km and the value of 
EDR is 3.3 x 10 47 erg /s, and the total energy deposited per second in the vicinity of the 
star (calculated by integrating EDR) is of the order of 10 49 erg. If we further increase the 
temperature to 20MeV the total energy deposited per second near the star is of the order 
of 10 51 — 10 52 erg, close to that of the energy liberated during GRBs. 

Let us try to understand the peculiar nature of initial increase in EDR with distance. 
From figure 3, for the isothermal disk the temperature To, seen by an observer at infinity, 
increases as the distance from the star increases. In this case T becomes hotter with distance 
and therefore it approaches T e fj (equal to Tc for isothermal disk) at large distance. The 
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EDR via neutrino pair annihilation is strongly dependent on disk temperature. Therefore 
let us calculate the pair annihilation of neutrinos emitted from the disk with TG. There 
is no definite theory governing the temperature gradient of the disk and we adopt the 
temperature gradient to be 1/r as the simplest acceptable model adopted by others for 
similar calculation [l6]. For comparison we also plot curves with TG = 1/r 1 / 2 . Due to 
variation of disk temperature, the temperature at the MPR (initial temperature from which 
our calculation starts) is much less than that of the central temperature. For Tc = lOMeV 
and TG = 1/r the T at MPR is less than 2MeV and the temperature falls very rapidly 
and reaches fractional values beyond 50km. The case of TG = 1/r 1 / 2 the situation is not so 
pessimistic. For the hotter star, T at the MPR is above 3MeV and does not falls to zero 
at large distances. 

Looking at figure 4, for the isothermal disk, we see that F(r) decreases with distance and 
falls to zero after 25km. For the disk with temperature gradient we see that F(r) initially 
increases and after about 8km it starts to fall, unlike the previous isothermal case where it 
decreases throughout. 

Initially for the isothermal disk the temperature dominates and EDR increases but after 
some distance (where F(r) decreases but temperature increases) F(r) becomes the domi- 
nating factor and the EDR decreases. The disk temperature observed at infinity is higher 
at outer region and the EDR increases with temperature, on the other hand the Doppler 
effect due to disk rotation reduces the EDR. The EDR calculated non relativistically at the 
MPR is half that of the relativistic case. 

Figure 5 shows the variation of EDR, for disk with TG, along the rotation axis of the 
star for two different central temperatures 5MeV and lOMeV. For Tc = lOMeV with 
TG = 1/r, the initial value of EDR is 10 m erg/s and it falls very sharply to zero at a 
distance of 50km. For Tc = 5MeV the fall is even sharper. For TG = 1/r 1 / 2 the fall in 
EDR is much slower and near the MPR the EDR, for central temperature Tc = lOMeV, is 
of the order of 10 A3 erg/s. For the colder disk the initial EDR is much smaller. These curves 
show that the EDR is very sensitive to temperature as well as to TG. At low r, temperature 
is the governing factor in determining the nature of EDR, therefore EDR decreases. As 
r increase F(r), which strongly decreases with r, becomes the dominant factor and EDR 
decreases further. This is why the EDR falls off much faster than the isothermal case. 

The slowly rotating star has a much smaller disk (smaller equatorial plane) and its MPR 
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is 4.7km, larger than the keplerian star. Figure 6, 7, 8, 9 depicts the various aspects of EDR 
for a slowly rotating star. The gross nature of the EDR for the isothermal case is more or 
less same, but the value of maximum EDR is slightly larger. The temperature plot is about 
the same, but the nature of F(r) is slightly different, the maximum value of F(r) being 
slightly larger. Due to this reason the EDR for the slow star is greater (twice) than that for 
the keplerian star. The rotational velocity seems to have some negative effect on the EDR 
and as the rotational velocity increases the EDR decreases. 

The value of EDR increases if we go inwards from the MPR. For any star (keplerian or 
slow) if we go close towards the centre, the EDR increases due to the increase in temperature 
of the disk. Near the centre the value of EDR for non relativistic calculation and that with 
GR effect along with TG becomes the same. But as we go outwards the value of the EDR 
for the GR calculation becomes less than the non-relativistic case. This is because, for GR 
calculation, the disk temperature falls as we go outwards along the disk. 

V. SUMMARY AND DISCUSSION 

In this article we have investigated the GR effect on the EDR of vv — > e + e~ reaction in 
a rotating neutron star described by CST metric. Here we have calculated the EDR for the 
neutrinos coming out of the equatorial disk and depositing energy along the rotation axis 
of the star, above the equatorial plane. The bending of the neutrino trajectories, and the 
redshift due to disk rotation along with gravitation has been taken into consideration. 

We find that for an isothermal disk, initially the EDR increases with distance and reaches 
a maximum value and then decreases with distance. This is due to the fact that there is a 
competition between the temperature and the Doppler effect due to disk rotation (charac- 
terized by F(r)). The temperature observed by the observer at infinity (T ) increases and 
at large distance becomes same as that of the effective temperature (T e ff) whereas F(r) 
decreases with increase in distance from the disk. In the case of disk with TG, the EDR 
falls off very quickly with distance. For TG — 1/r the slope with which EDR falls is much 
greater than that for TG = 1/r 1 / 2 . The F(r) increases initially but it also decreases after 
some distance and falls very near to zero at about 50km. Both the isothermal and TG 
picture shows that the EDR is a very sensitive function of temperature, and as temperature 
decreases the EDR also decreases. The slow star shows more or less the same nature only 
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with the fact that EDR for the slow star is about twice that for the keplerian star. The 
rotational effect reduces the EDR. 

To summarize, we find that the EDR is very sensitive both to the temperature and TG 
of the disk. The deposition energy is contributed mainly by the neutrinos arriving from the 
central region where the temperature is higher. The maximum energy is deposited near the 
surface. As we move outwards the deposition energy reduces and is smaller than that of the 
non relativistic calculation. 

The total energy deposition from a neutron star would also depend on the off-axis con- 
tributions. So to get an estimate of the total deposition energy, one needs to do a detailed 
numerical simulation taking into account all possible contributions. Such a calculation is in 
progress and we hope to report it in the near future. 
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FIG. 2: Variation of EDR with distance along the rotation axis for the keplerian star having 
isothermal disk with three different central temperature. 
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FIG. 3: Variation of temperature To with distance along the rotation axis of the keplerian star. To 
is plotted for both isothermal disk and for disk having two different temperature gradient. 



12 




FIG. 4: Variation of F(r) with distance along the rotation axis of the keplerian star. F(r) is 
plotted both for isothermal disk and disk with temperature gradient. 
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FIG. 5: Variation of EDR with distance along the rotation axis for the keplerian star. The plot 
have been done with disk having two different TG and the curves are plotted with two different 
central temperature. 
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FIG. 6: Variation of EDR with distance along the rotation axis for the slow star having isothermal 
disk with three different central temperature. 
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FIG. 7: Variation of temperature To with distance along the rotation axis of the slow star. To is 
plotted for both isothermal disk and for disk having two different temperature gradient. 
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FIG. 8: Variation of F(r) with distance along the rotation axis of the slow star. F(r) is plotted 
both for isothermal disk and disk with temperature gradient. 
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FIG. 9: Variation of EDR with distance along the rotation axis for the slow star. The plot have 
been done with the disk having two different TG and the curves are plotted with two different 
central temperature. 
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